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Abstract 

We study the stability properties of periodic solutions to the Nonlin- 
ear Schrodinger (NLS) equation with a periodic potential. We exploit the 
symmetries of the problem, in particular the Hamiltonian structure and 
the U(l) symmetry. We develop a simple sufficient criterion that guaran- 
tees the existence of a modulational instability spectrum along the imagi- 
nary axis. In the case of small amplitude solutions that bifurcate from the 
band edges of the linear problem this criterion becomes especially simple. 
We find that the small amplitude solutions corresponding to the band 
edges alternate stability, with the first band edge being modulationally 
unstable in the focusing case, the second band edge being modulation- 
ally unstable in the defocusing case, and so on. This small amplitude 
result has a nice physical interpretation in terms of the effective mass of 
a particle in the periodic potential. We also consider, in somewhat less 
detail, some sideband instabilities in the small amplitude limit. We find 
that, depending on the Krein signature of the collision, these can be of 
one of two types. Finally we illustrate this with some exact calculations 
in the case where V(x) is an elliptic function, where all of the relevant 
calculations can be done explicitly. 
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1 Introduction 



In this paper we consider the stability of standing wave solutions to the NLS 
equation with an external potential: 

i^ t = -^ xx ±e\^ + V(x)iP, (1) 

where the potential V(x) has period 1. We assume that we have a standing wave 
solution tp(x, t, e) = exp(— iu>(e)t)<f)(x, e), where <j)(x, e) is real valued and either 
periodic (cf>(x + 1, e) = <f>(x, e)) or antiperiodic {<j>(x + 1, e) = —<f>(x, e)). For cer- 
tain perturbative results we will also assume that the standing wave bifurcates 
continuously from the Floquet-Bloch eigenfunctions in the usual way. In other 
words we assume that <f>(x, 0) is an eigenfunction of the periodic Schrodinger 
operator 

fMt> = u{0)<f>=~<f> xx + V(x)<f> (2) 

corresponding to a band edge, and <p(x, e) e C 2 ([0, 1] x (—5, 5)) 

We find that the combination of the Hamiltonian structure and the U(l) 
symmetry dramatically simplifies the structure and possible bifurcations of the 
spectrum of the linearized operator. We find a simple condition on the L + op- 
erator which guarantees the existence of a modulational instability. In the case 
of weakly nonlinear standing waves, when e is small, we apply a perturbation 
argument to show that the lower band edges are modulationally unstable in the 
focusing case, while the upper band edges are modulationally unstable in the 
defocusing case. 

We also study, in somewhat less detail, some bifurcations which occur for 
non-zero values of \x when a certain eigenvalue degeneracy condition is met. We 
find that this bifurcation generically leads to the emergence of complex eigenval- 
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ues, and thus instability. Depending on the Krein signature of the unperturbed 
eigenvalues this may or may not lead to the opening of a gap along the real axis. 

2 Fundamentals of Systems with Symplectic Struc- 
ture and Periodic Coefficients 

In this section we review some results of Floquet theory for the special case in 
which the equations admit a Hamiltonian formulation. For more details see the 
text of Yakubovich and Starzhinskii|18|. 

Throughout this paper we will apply the following notation. U(x) will denote 
the solution operator of the periodic ODE governing stability. The matrix M 
will denote the monodromy, or period, map U(l). We will have occasion to 
consider the spectral properties of the second order operators L± individually. 
The matrix m will denote the mondromy matrix of the second order problem 
associated to L± . For each of these quantities the dependence on the parameters 
fi, e will be generally be suppressed unless it is necessary for clarity. Similarly 
K±(fi) will denote the Floquet discriminants of the full stability problem, while 
k(fi) will denote the Floquet discriminant of the second order problem associated 
with L±. 

2.1 The Lyapunov-Poincare theorem 

To begin we assume a Floquet problem of the following form: 

U x =JH(x)U U(0) - I 2JV x2JV H(a; + l)=H(x) H* = H (3) 

where, for simplicity, we have set the period to one. Here l27Vx27V is the 2Nx2N 
identity matrix and J is a skew-symmetric matrix. We will assume that J 2 = 
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— I since the above problem with a non-singular J can always be mapped to 
one of that form. The monodromy matrix M is defined to be the period map 
M = U(l). It is easy to see that U(x), and thus M, satisfies the relation 

U'JU=J, (4) 

so that U(a:) is a symplectic 3 matrix. From this it follows that det( U) = 1, and 
more generally that the characteristic polynomial P(X) = det( U — A I) satisfies 

det(U-AI) = dct( U* — AI) = dct( J( U f — A I) J*) (5) 
= dct(U- 1 -AI) = A 2iV dct(U- 1 )det(U-A- 1 I) (6) 
= A 2A 'P(A- 1 ) (7) 

thus the polynomial is palindromic - if P(X) = J2 a 3 X then 
This result is known as the Lyapunov-Poincare theorem JS] . 

This symmetry implies that the problem of finding the roots of the 2N th 
degree polynomial P(X) can be reduced to finding the roots of an N th degree 
polynomial P{z) by means of the transformation z = A + A . For the case of 
the stability problem for standing wave solutions to the nonlinear Schrodinger 
equation the monodromy matrix is 4 x 4 (N = 2) In the case the characteristic 
polynomial of the monodromy matrix takes the form 

P(X) = l + aX + bX 2 + aX 3 + A 4 = 0. (8) 

Under the conformal map the roots of this polynomial are mapped to the roots 
of the second degree polynomial 

P(z) = z 2 + az+ {b- 2) = (9) 

3 Sometimes called J-unitary 



4 



with A = Z± V^ 2 4 ; an d the characteristic polynomial admits an explicit factor- 



ization into two second degree polynomials 



P(A) = (1 - K+X + A 2 )(l - K-\ + A 2 ) 



(10) 



where the K± are the following algebraic functions of a, b: 



K± = 



-a ± Va 2 - 46 + 8 



(11) 
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The functions if± will be refered to as the Floquet discriminants. Since the 
conformal mapping z = A + j takes the unit circle to the real interval [—2,2] it 
follows that the monodromy matrix has two eigenvalues (counted by algebraic 
multiplicity) on the unit circle if K + lies in the real interval [—2,2], and two 
more if K_ e [-2,2]. 

It is convenient to express the K± in terms of the invariants of M. Since 
a = - Tr( M), b = -\ Tr( M 2 ) + \ Tr( M) 2 the Floquet discriminants arc given 



This leads to our first observation 

Lemma 1 The possible (algebraic) multiplicities of the eigenvalues of a 4 x 4 
monodromy matrix are (1, 1, 1, 1) , (1, 1, 2), (2, 2) and 4. The conditions on the 
monodromy matrix M which produce eigenvalues of higher multiplicity are as 
follows: 



by 




(12) 
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Condition on K± 


Condition on M 


Root Multiplicities 


K+ = ±2 or K- = ±2 
Simple Band Edge 


Tr(M 2 ) = (Tr(M)±2)) 2 


(1,1,2) 


K+ = ±2 and K- = t2 
Repeated Band Edge 


Tr(M) = 
Tr(M 2 ) = 4 


(2,2) 


K+ = K_ ^ ±2 
'Accidental' Degeneracy 


2 Tr(M 2 ) = (Tr(Af)) 2 -8 
Tr(M) ^ ±4 


(2,2) 


K+=K-= ±2 
Double Band Edge 


Tr(M) = ±4 
Tr(M 2 ) = 4 
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Proof: A straightforward calculation from the explicit factorization given in 



It is "worthwhile to consider the codimensions of the above possibilities. Since 
we have a one parameter family of mondromy matrices, parameterized by the 
spectral parameter /j,, we should generically expect to see possibilities 1 and 3, 
which require only one condition, but not possibilities 2 and 4, which require 
two independent conditions be satisfied. We shall see that possibility 4 always 
happens at fj, = for small amplitude waves, and is forced by symmetries. It is 
this degeneracy which leads to the modulational instability. It is this situation 
that we will consider in most detail. We will also consider possibility 3, in 
somewhat less detail. We will not consider possibility 2, since it is non-generic, 
nor possibility 1, since it is not very interesting - one expects that the band 
edges will typically move under perturbation. 
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2.2 Krein Signature 

There is an extensive theory of the structural stability of symplectic matrices due 
to Krein and collaborators. We shall give only a brief summary of this theory 
here. It is clear that, if a symplectic matrix M has k distinct eigenvalues on 
the unit circle it follows from continuity that a nearby symplectic matrix M 
must also have k distinct eigenvalues on the unit circle. Thus eigenvalues of 
symplectic matrices can only leave the unit circle via collisions. What is less 
clear is that only certain collisions can lead to pairs of eigenvalues leaving the 
unit circle. The quantity which distinguishes such "dangerous" collisions is that 
of the Krein signature. The Krein sign r\ of an eigenvector v of a symplectic 
matrix M is defined to be 

77 = sgn(<v 3v>), (13) 

while the signature (p, q) of an r— dimensional eigenspace is defined as follows: p 
(resp. q) is the number of linearly independent eigenvectors with positive (resp. 
negative) Krein sign. 

The fundamental stability result for symplectic matrices says that a sym- 
plectic matrix is structurally stable (i.e. the number of eigenvalues on the unit 
circle does not change under perturbation 4 ) if all eigenspaces corresponding to 
eigenvalues on the unit circle have a definite Krein signature. Further this result 
is tight: if an eigenspace is of indefinite Krein signature then a generic pertur- 
bation will cause roots to move off of the unit circle. The first result is known as 
Krein's theorem, while the second is known as the Krein-Gelfand-Lidskii strong 
stability theorem. For precise statements, as well as proof, see Yakubovich and 
Starzhinskii 18 . 

4 Here perturbations are always assumed to preserve the symplectic nature of the matrix. 

7 



2.3 Hill's Equation 



In this section we state some properties of the Hill's equation 



tpxx + Q(x)lp = flip 



(14) 



that will be useful in the sequel. We will also define a quantity that will be 
useful in the analysis to follow. 

In the remainder of this section m is the (2 x 2) monodromy matrix associated 
to the Hill's equation l(TIjl: 



We begin with some well-known facts. The spectrum of the above eigen- 
value problem consists of the union of the set of intervals spec(-ff) = [/to, Mi'] U 
[/^'j Mi] U [a*2j M3 7 ] ■ ■ ■ j where the fa are the roots of k(fi) = 2 and the fa' are 
the roots of fc(/i) = —2. The points {fa} are the periodic eigenvalues, while the 
points {fa'} are the anti-periodic eigenvalues. Together these comprise the band 
edges. It is always true that fan < fai+i' and fjf 2i < fai-i, so the bands are 
always nontrivial. It can happen that faii-i — A*2i or A*2i-l — A*2j corresponding 
to a closed gap. This is commonly referred to as a double point. In the interior 
of each band there are two quasi-periodic solutions, while at the band edge there 
is one periodic or antiperiodic solution and one solution which grows linearly 
unless the band edge is a double point, in which case there are two periodic or 
antiperiodic solutions. The derivative of the Floquet discriminant is non-zero in 



rn = 




(15) 



with k(fa) = Tr(m), and j is the standard skew-symmetric form 



J 




(16) 
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the bands, has exactly one zero in each gap, and is nonzero on the band edges 
unless the band edge is a double point, in which case the derivative of Floquet 
discriminant has a simple zero. Proofs of these facts can be found in Magnus 
and Winkler [TT]. 

In the first lemma we introduce a quantity which will play an important 
role in the perturbation analysis. The sign of this quantity will be important in 
determining the modulational stability of the standing wave. 

Lemma 2 Define the quantity <r(/x) as follows: 

a(fi) = Tr(jm(//)) (17) 

Then the following hold 

• o~(/i) is non-zero on the interior of each band. 

• vanishes at a band edge iff the band edge is a double point. 

• sgnH/j)) = -sgn(fc'Ou)) = -sgn(Tr(fs)) if ^ e spe c(H). 

Proof: On the interior of a band it holds fc 2 (/i) — 4 < 0. If we set r/i = 
m = ^2(1, m)> Vi = an d i 2 = ^(i)/"), </l e«. 

fc 2 ( M )-4= (tji-^) 2 +4^»h, (18) 

which implies that r/^r/2 < 0. 5j/ iis definition, a(fi) = — 77^ +772 = — sgn(7]' 1 )(|?7' 1 1 + 
I ] ) - implies that rj[ 7^ and 772 7^ 0. TTius it follows that c(/i) is 

nonzero in the interior of each band. 

Let a band edge be a double point. It follows that k 2 (p) — 4 = and moreover 
k'(p) = 0. Following 11 If one can show that 

fc'(/x) = sgnfoixKi / V 2 + \ m \ I i>! ± 2 X M^| / V'i^), (19) 



where ipi, -02 are the two linearly independent solutions of 

ipxx + (M - Q(aOW = °> ( 2 °) 

that satisfy the initial conditions ipi(0,/i) = 1, ^(0,/it) = 0, ^2(0, /*) = 0, and 
^2(0, /x) = 1. FFe Mserf again fltf) and, in particular, the inequality r\' x r\2 < 0. 
By Cauchy- Schwartz it follows that sgn(?7j) = 0, thus rj-y = 0. But then, also 
rji — r(i — 0, which in turn forces r\i = 0. 

The other implication is proved as follows: o~(/z) — implies = fc 2 (/i) — 4 = 
(7/1 — 77 2 ) 2 + 4?72, 50 772 = ??i = and 7/1 — r/ 2 =0. But then clearly k'(p) = 0. 

Finally, again in it was shown that if rf^ 7^ then k' (n) and 77^ have the 
same sign. The result follows, since from the first part of this lemma, o~(p) and 
r{ x have opposite signs. 

This quantity is related to the Krein sign, and actually agrees with the 
Krein sign of the eigenvalue of the monodromy matrix in the upper half plane. 
However the main utility of this quantity is that it allows one to compute the 
sign of the off-diagonal piece of the Jordan normal form at a band edge, as in 
the lemma below. 

Lemma 3 At a band edge the monodromy matrix m has the following Jordan 
normal form: 

( ±1 trfa 

m = o 

. ±1 

where o is a proper orthogonal matrix: oo* = I, det(o) = +1 

Proof: It is known that at a band edge that is not a double point ±1 is 
a7i eigenvalue of algebraic multiplicity two, and geometric multiplicity 1. The 
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Jordan normal form implies that 



±1 K 



\ 



m = c/ 



o = o mo 



(22) 



V 



o 



±i 



J 



where o is orthogonal and can be chosen to have determinant +1. It remains 
to be checked that K = cr(/i). One easily observes that'} is a rotation and thus 
commutes with o. This implies that a = Tr(jm) = Ti^ojc^m) = K. If o is 
chosen to be on the other connected component of the orthogonal group j and o 
anti-commute, and the sign of the off-diagonal element is reversed. 

The value of cr at a band edge, in particular the sign, is important in de- 
termining the stability of small amplitude standing wave solutions. Note that 
the Krein signature of an eigenvalue of the monodromy matrix in the upper 
half-plane is equal to the sign of a, and thus has the opposite sign from the 
derivative of the Floquet discriminant. 

3 Modulational Instability of Standing Waves 
3.1 General Results 

We assume that there exists a C 2 family of standing wave solutions ip = 
e~ w ^ t 4> s tand{x,e) to the NLS equation with a periodic potential 



which is either periodic 4>{x+ 1) = <f>(x) or anti-periodic <f>(x+ 1) = —<p(x). The 
eigenvalue problem governing the linearized stability is given by 




L-5 = 
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where the operators L± are given by 

L- = ~d xx + V(x) + e\(j) s tand\ 2 {x,£) - w(e) 
L+ = -^d xx + V(x)+3e\(f> stand \ 2 (x,e)-uj(e) 

The above eigenvalue problem has a Hamiltonian formulation for real /x. The 
Lyapunov-Poincare theorem of the previous section implies that the characteris- 
tic polynomial of the monodromy matrix is palindromic, P(X) = A 4 P(A _1 ), for /j, 
on the real axis. Also note that for arbitrary values of e the spectrum is symmet- 
ric about the real and imaginary axes, since the eigenvalue problem is invariant 
under the transformations /j — > — fx,p — ► — p, g — > g and /i — > /2,p — > p, g — > g. 

We begin by noting a few simple properties of this eigenvalue problem, and 
the associated monodromy matrix: 

Proposition 1 The Floquet stability problem for an NLS standing wave has 
the following properties: 

• 11 = is a periodic (anti-periodic) eigenvalue. 

• M(/i) is an entire function of [i, of order |. 

• The monodromy matrix M is symplectic for all fi e C. 

• The Floquet discriminants fc±(/x) are analytic functions of /i away from 
the branch points where 2Tr(M 2 ) - Tr( M) 2 = 8. 

• At ji = the characteristic polynomial of the monodromy matrix has the 
following special form: 

P(A)| M=0 = 1- Tr(M)A + (2Tr(M)-2)A 2 - Tr( M)A 3 + A 4 (24) 
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Proof: 



That /i = is always an eigenvalue follows from Noether's theorem and the 
phase invariance of NLS. The corresponding eigenvector is p = 4> s tand{x), q = 0. 
The fact that M is entire follows from standard arguments - see for example 
the text of Sibuva[T7|. Note that 2 Tr( M 2 ) - Tr(M) 2 = 8 is also an entire 
function of fractional order, and thus must have a countable number of zeros. 
The fact that M is symplectic for real fj, follows from the results cited in the 
previous section. To see that this in fact holds for all fi E C simply note that 
J — M* J M is an entire function that is zero on the real axis, and thus must be 
identically zero. A simple division shows that the polynomial A — 1 divides the 
polynomial A 4 + a A 3 + 6A 2 + a A + 1 if and only if b + 2a + 2 = 0, which proves 
the last part. Note that this same calculation shows that if A = 1 is a root it is 
necessarily of multiplicity 2 or 4. 

Lemma 4 If is not a periodic eigenvalue of the L + operator the Floquet dis- 
criminants K±(fj.) are analytic in a neighborhood of /J, = 0. 

Proof: For /i = the stability problem decouples, and the monodromy matrix 
takes the block diagonal form 



where m± are the monodromy matrices associated with L±. It follows from 
the previous lemma that A = 1 is an eigenvalue of m with multiplicity 2. A 
short calculation using the fact that 2x2 matrices satisfy Tr(m 2 ) = Tr(m) 2 — 
2det(m) shows that 2 Tr( M 2 ) - Tr( M) 2 + 8 = (Tr(m_) - 2) 2 . If is not a 
periodic eigenvalue of L + then 2 Tr( M 2 ) - Tr( M) 2 + 8 ^ and K ± arc analytic 
in a neighborhood of fi = 0. 




(25) 
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Theorem 1 For general e a sufficient condition for the existence of a modula- 
tional instability spectrum is that fi = is in the interior of a band of the L + 
operator. 

Proof: From the previous lemma if yu = is in the interior of a band of the 
L + then the Floquet discriminants K ± (/j.) are analytic in a neighborhood of the 
origin. From the fact that the coefficients of the characteristic polynomial are 
invariant under the transformation fj, — > — /j, it follows that K±(/j,) are even, and 
thus are real on an segment of the imaginary axis containing the origin. Since 
if-(O) G (—2,2) it follows that K_(jj) is real and G (—2,2) on some segment 
of the imaginary axis containing the origin, and thus there is a modulational 
instability. 

Remark: The same argument guarantees the existence of a modulational 
instability if cither K(0) = +2 and K"(0) > or K(0) = -2 and K"(0) < 0. 
The calculation of the sign of the second derivative appears to require a difficult 
second order perturbation calculation. 

3.2 Perturbative Results for weak nonlinearity 

In this section we study the Floquet spectrum for small amplitude standing 
waves of the nonlinear Schrodinger equation with a periodic potential, with a 
particular emphasis on the behavior near \i = 0. We shall see that the U(l) 
phase invariance of the NLS forces a four-fold degeneracy of the eigenvalues 
of the unperturbed monodromy matrix at fi = 0. Under perturbation this 
degeneracy can lead to the birth of a spine of continuous spectrum lying along 
the imaginary axis. Whether or not such a spine is born is determined by 
the relative sign of the nonlinearity and the quantity a defined in the previous 
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section, and the length of the spine is determined by the magnitude of a. 
When e = the operators L± are equal and are given by 

L-{0)=L+{0) = ~d xx + V(x)-u(0) (26) 

and the resulting eigenvalue problem is self-adjoint. It is easy to see that in this 
case the spectrum of the stability problem is given by spec(L_) U spcc(— L_). 
In this case it is also straightforward to calculate that the monodromy matrix 
M of the full stability problem takes the block diagonal form 

( m (M) \ 
M = (27) 

where m(±/x) is the monodromy matrix for the problem 

L_(0)^> = ±juV- ( 28 ) 

From this block diagonal form it follows that for e = the invariants of the 
full monodromy matrix can be expressed in terms of the Floquet discriminants 
fc(±/x) of the second order problem via 

Tr(M( M )) = Tr(m(/i)) + Tr(m(-/i)) = k(fi) + fc(-/i) 
Tr(M 2 ( M )) = Tr(m 2 ( M )) + Tr(m 2 (- M )) = k 2 (p) + k 2 {- - A. 

Here we have used the fact that 2x2 matrices satisfy the identity Tr(m 2 ) = 
Tr(m) 2 — 2det(m). From this it follows that the Floquet discriminants K±(fi) 
for the full problem can be written in terms of the Floquet discriminant of L 
via 



Obviously this could be simplified to K±(p) = fc(±/x) however we do not do 
this, since it obscures the degeneracy at the points where = /i). 
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The U(l) symmetry of the NLS equation implies that \i = is a band-edge 
eigenvalue of L . Thus for e = 0, fi = the monodromy matrix has has a 
single eigenvalue ±1 of multiplicity four. This corresponds to the last entry 
in the table in figure 1. This seemingly non-generic bifurcation is forced by 
the phase-invariance symmetry along with the symplcctic nature of the mon- 
odromy matrix, and can give rise to the modulational instability at non-zero 
wave amplitudes. 

In the next lemma we present a normal form calculation for the Floquet 
discriminants in a neighborhood of /i = 0, e = 0. The calculation is particu- 
larly simple due to some additional symmetries, which dramatically reduce the 
number of coefficients which need to be computed. We present the calculation 
for periodic band edges: the calculation for the anti-periodic bands edges is 
identical except for some sign changes. In this latter case we merely state the 
final result. 

Lemma 5 The Floquet discriminants of the monodromy matrix have the fol- 
lowing normal form in a neighborhood of e = 0, fi = 0: 

K ± {n) = 2 + + 2eo- <0f> +Ex ± y / (fc'(0) A i) 2 + (2ae«j>\>f + E 2 

Ei -o(e, M 2 ) 
£ 2 =o(e 2 , M 2 ) 

where k(fi) is the Floquet discriminant for the L_(0) operator. 

Proof: This perturbation calculation is somewhat tedious but can be made 
simpler by the use of some of the previously derived identities. From the last 
part of Lemma 4 we can express the e derivatives of the coefficients of the 
characteristic polynomial at the origin in terms of Tr(M|/i = 0). To compute 
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the jU derivatives of the coefficients of the characteristic polynomial at the origin 
we use the fact that a(/z, 0) = — (fc(/i) + k(— fi)), 0) = k(fi)k(— /j,) + 2. The 
mixed partial vanishes at the origin since the o, 6 coefficients are even 

functions of /x for all e. Note that, because of the square root branch point, it 
is necessary to carry the expansion out to second order to get what amounts to 
a first order result - the local normal form is a cone. 

We present the calculation of d ^- T( - ^ \^-o) g rst q^e problem L_?/j = 
has a periodic solution that we denote by ipi(= ^ 'stand), and a linearly growing 
solution that we denote by ip2- We choose ^1,2 to form an orthogonal right- 
handed coordinate system: 

V>i(0) = cos(6>) = sin(0) 

^ 2 (0) = - sin(6>) V 2 (0) = cos(6>) 

Since \i = is a band edge for the L operator we have the following expressions 
for the period map: 

$j(l) = ^j(0)+ o^i (0) = CTCOs(0)-sin(6») 
^a(l) = crsin(6») + cos(0). 

Note that a is non-zero as long as the band edge is not a double point. It 
is convenient to define a second set of solutions ^>i ; 2 which satisfy L_0 = 
along with the initial conditions 0i(O) = 1,^(0) = and </>2(0) = 0, </>' 2 (0) = 
1. These are obviously related to the ^1,2 by <p\ = cos(0)V>i — sin(0)^>2 and 
<p2 — sin(#)^>i + cos(6)ip2- These functions form a natural basis in which to do 
perturbation theory on the L + operator. As is usual in Floquet theory we must 
construct a fundamental set of solutions 4>i , 4*2 to 

L+01,2 =L_0 l!2 + 2e|V| 2 0i,2 = (30) 
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that satisfy the boundary conditions 



0i (0) 



1 #(0) = 



02(0) 



0' 2 (O) = 1 



A straightforward perturbation calculation gives the following expressions for 
the fundamental set of solutions to L + = : 



Using the fact that "02(1) = "02(0) + aipi{0) gives the following expression for 
the trace of the monodromy matrix to the leading order in e : 



0i(l) = 1 (l)-4e(cos 2 (f?)(<^ 2 >-a<^>) 

+ sin 9 cos 6(a < ^f<0 2 > + < i>\ > - < i>li>l >) 

- sin 2 (0) < iftfo >) 

2 (1) = 2 (l)-4e(~cos 2 (0) <^ 2 > 

+ sin(0) cos(0) < i>\i>\ > -a < ^f<0 2 > - < i>\ >) 

+ sin 2 9(< ipfip 2 > -o- < Vi >) 



where < / >= J Q f(y)dy. From this it follows that the monodromy matrix in 
leading order of e is given by 



Tr(M| M = 0) = 1 (1)+0 2 (1)+0 1 (1)+0 2 (1) = 2(0 1 (l)+0' 2 (l))+4e ( 7 <tf> +0(e 2 ) 

(31) 



It is easy to calculate that a(n, 0) = — 4 — fc"(0)/i 2 + o(/i 2 ), and that 0) = 
6 + (2£;"(0)- (fc'(0)) 2 )^ 2 + o(^ 2 ). This, together with the vanishing of the mixed 
partial and the explicit formula gives the above result. 
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Remark From the proof of the preceeding lemma, and the use of the time- 
invariant Hamiltonian energy functional 

n = J Ql^l 2 ± f - V(xM 2 ^j dx, (32) 



it follows that 

dTrM 



(0,0) ut 



(33) 

s=0 



ck 

Thus, the instability criterion can be expressed in terms of the signs of a and of 
the derivative of the energy functional. This is a feature that appears regularly in 
the study of stability of nonlinear waves: that the sign of the derivative of some 
conserved quantity is a proxy for the index of some linearized operator. For some 
results of a similar flavor in a variety of different contexts (both conservative 
and dissipative) see [T3l 151 171H21 15| . 

Theorem 2 For e small and positive (focusing NLS) and solutions bifurcating 
from a lower band edge that is not a double point, or for e small and negative 
and solutions bifurcating from an upper band edge that is not a double point 
there exists a band of spectrum along the imaginary axis. 

Proof: The proof follows from the preceeding lemma and the same ana- 
lyticity argument as the first theorem. From the explicit factorization of the 
characteristic polynomial it follows that, for e fixed, the Floquet discriminants 
K± (ji) are analytic functions in the cut plane with branch points at the points 
where 

2Tr(M 2 ) — Tr(M) 2 + 8 = (34) 
From the preceeding Lemma it follows that at /i = this quantity reduces to 

(2Tr(M 2 ) - Tr(M) 2 + 8) | M = = (Tr(M)| Al=0 — 4) 2 = 16a 2 e 2 <^> 2 +o(e 2 ) 

(35) 
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and thus this quantity is nonzero for small but nonzero e: the fact that the band 
edge is not a double point guarantees that a ^ 0. It follows that the Floquet 
discriminants K± are analytic in a neighborhood of \x = and satisfy 

K+(0) = 2 + Aev < i>\ > +0(e 2 ) 
K-(0) = 2. 

As in the first theorem the fact that K + (p) is an even function implies that 
(i = is a critical point, and K + is real on a segment of the imaginary axis 
containing the origin. Recall that for odd numbered bands the periodic eigen- 
value is the lower band edge, and a is positive, and for even numbered bands 
the periodic eigenvalue is the upper band edge, and a is negative. For e < 
(defocusing) and odd numbered bands, or for e > and even numbered bands 
the perturbation result implies that K + (0) £ (—2, 2), so there exists an interval 
on the imaginary axis on which K + (fi) is real and £ (—2,2). The analogous 
calculation for solutions bifurcating from the anti-periodic eigenvalues shows 
that for e < (defocusing) and even numbered bands, for e > and odd num- 
bered bands K + (0) £ (—2,2), and there exists a spine of spectrum along the 
imaginary axis. Thus, for e small and defocusing nonlincarity standing wave so- 
lutions which bifurcate from the lower band edge are modulationally unstable, 
while for focusing nonlinearity standing wave solutions which bifurcate from the 
upper band edge are modulationally unstable. Said more simply the solutions 
which bifurcate from the lower band edges are modulationally unstable in the 
defocusing case, and the solutions which bifurcate from the upper band edges 
are modulationally unstable in the focusing case. 

Note that if K±(Q) = there exist additional arcs of spectrum emerging 
from the origin into the complex plane. For e sufficiently small the local normal 
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form calculation guarantees that the second derivative is non-vanishing, and 
this does not occur. It is a possibility for larger e, however. 

4 Sideband Instabilities 

In this section we consider the possibility of sideband type instabilities. This 
case is somewhat more difficult than the case of modulational instabilities, since 
there are fewer symmetries, and our results are somewhat less detailed. In 
addition to the analyticity arguments of the previous section a major tool will 
be the Krein signature. 

First we note the following lemma 

Lemma 6 The spectrum of the stability problem consists of a union of contin- 
uous curves in the complex plane. The possible endpoints of the curves are band 
edges K±(fi) = ±2, branch points a 2 (/z) — 46(/i) + 8 = (or K + (p) = K^(n)) 
or critical points K'±([i) = 0. 

Proof: As before K±(/j.) are analytic away from the (isolated) points where 
a 2 (/i) — 46(/i) + 8 = 0. Suppose \i is a point in the spectrum that is not a 
branch point, band edge or critical point. We assume for the sake of argument 
that K + (ji) e (—2,2). Since K + is analytic in a neighborhood of ^ and fi is 
not a critical point of K + an appeal to the implicit function theorem shows the 
existence of a unique curve through fi along which Xm(K + ) = 0. Since \i is not 
a band edge continuity implies that lZe(K + ) G (—2, 2) in some neighborhood of 
11. 

We now consider the possibility of sideband instabilities arising from points 
where K + ([i) = K_([i). First note that in the small amplitude limit, when e = 0, 
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we have a 2 — 46+8 = (k(pi) — k(— pi)) 2 so the branch points are degenerate. Under 
perturbation this degeneracy is expected to break, and generically it can do so in 
one of two ways, which are generally referred to as the avoided collision and the 
open gap. These possibilities are illustrated in figure |SJ It is rather tedious to 
compute the normal form for this case, since one lacks the even symmetry, but 
consideration of the Krein-Gelfand-Lidskii theorem can reduce the possibilities. 
In the case of a collision of two eigenvalues with the same Krein sign the open 
gap is forbidden, since that corresponds to eigenvalues leaving the unit circle. 
In the case of two eigenvalues with the opposite Krein sign one expects that 
under perturbation the eigenvalues should leave the unit circle. We formalize 
this observation as a theorem. 

Theorem 3 In the neighborhood of an accidental degeneracy (pi — pi*,e = 0) 
the Floquet discriminants have the following normal form 

K ± (pi) = k(pi*) + (k'(pi*) - k'(-pi*))(pi - pi*) +ae + E x ± 
(k'(pi*) + k'(-pi*)f (pi - pi*) 2 + 7 e* + E 2 
E 1 ^0(e 2 ,(pi~pi*) 2 ,e(pi-pi*)) 
Ex = 0(e 3 , (pi - pi*) 3 ,e(pi - pi*) 2 ,e 2 (pi - pi*)) 

Further, ifsga(k'(pi*)) = sgn(fc'(— pi*)) the quadratic form (k 1 (pi*)+k' (—pi*)) 2 (pi— 
pi*) 2 + (3e(pi — pi*) + 7e 2 is nonnegative. 

A straightforward, though slightly tedious, perturbation argument shows that 
■§l(a 2 — 46 + 8) j e =o = 0, so that this quantity is locally quadratic, as well as 
having a vanishing mixed partial. The Krein signs of the unperturbed eigen- 
values are the same as the signs of k'(pi) and k'(-pi) respectively. In the case 
where these signs are the same the Krein-Gelfand-Lidskii theorem guarantees 
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that, under small perturbations, the eigenvalues of M remain on the unit circle. 
It is clear that a necessary condition for the eigenvalues to remain on the unit 
circle is non-negativity of the quadratic form, thus the sign condition guaran- 
tees non-negativity of this form. Note that this theorem does not preclude the 
following possibilities: in the case of like Krein sign collisions it is possible that 
the quadratic form has a zero eigenvalue, or that the intersection persists under 
perturbation. In the case of collisions of opposite Krein sign a naive applica- 
tion of the Krein-Gelfand-Lidskii theorem gives no information. We conjecture 
that, in this case, the signature of the quadratic form is equal to the signature 
of the subspace. This would imply that that intersections of opposite Krein 
signs always open to a gap, while intersections of the same Krein sign always 
open to an avoided collision. Very preliminary numerical evidence has supports 
this conjecture, but unfortunately it appears that the only way to verify this 
analytically is via a very tedious perturbation argument. 
An obvious corollary of this is the following: 

Corollary 1 If sgn(fc'(/i*)) = sgn(fc'(— fi*)) and 7 > 0, or if sgn(fc'(/i*)) = 
— sgn(fc'(— /j*)) and 7 < then there exists a band of spectrum off of the real 
axis. 

From the local normal form it is an easy calculation that, in either case, 
for sufficiently small e the Floquet discriminants K± must have a critical point 
in the neighborhood of fi* , and that a 2 — 46 + 8 7^ at this critical point. 
From the analyticity arguments of the previous section this guarantees that K ± 
are real and G (—2, 2) in some neighborhood of the critical point. Again one 
expects that in the case of a like Krein sign collision one should generically have 
7 > 0, leading to the avoided collision and a loop of spectrum opening into 
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the complex plane, but there seems to way to show this in any particular case 
without actually calculating 7. 



5 Explicit Examples and Numerics 

In this section we present some examples. We will primarily be working with 
known exact elliptic function solutions. We consider the nonlinear Schrodinger 
equation 

iA = -^ xx + V sn 2 (x,k)^±\^\ 2 ^ (36) 
This equation has a one-parameter family of exact solutions given by 

ip(x,t) = r(x)e- tut+te{x) (37) 
r 2 {x) = Asn 2 (x,k) + B (38) 

- *f ^7) (39) 
A = -(Vo + fc 2 ) (40) 
c 2 = *(^^) + (41) 

Where B e (-00, -fc 2 ) U (^T^,Vo + k 2 ^ for the focusing sign and B e 
(Vo + fc 2 ), — V °^i k ^ U (— fc 2 , 00) for the defocusing sign. These solutions rep- 
resent nonlinear stationary states which bifurcate from the linear Bloch states. 
We are primarily interested in the solutions which bifurcate from the band 
edges. These correspond to the boundaries of the above regions of validity. In 
the focusing case we have 



ro(,M V ^+ fc2 dn(x,fc) u = -1 - ^ + ^ \ - + > „ , 4 2) 



n{x) = VVo + k 2 cn(x, k) uj=^-V -k 2 v + k 2 > (43) 
r 2 (a:) = V-(Vb + fc 2 )sn(a;,fc) w = V + fc 2 < 0, (44) 
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while in the defocusing case we have the analogous solutions 



ro(x) 




k 2 Vr, k 2 

-dn(x,k) u = -l-j°+!L V Q + k 2 < 







(45) 



n{x) 




(46) 



T2{x) 




(47) 



The solutions all bifurcate from the the linear Bloch states at Vq + k 2 = 0, with 
the dn, cn solutions existing on one side of the bifurcation and the sn solution on 
the other side. It is easy to check that the spectrum of the L_ operator in each 
of these cases is the 1-gap Lame operator plus some constant which differs in 
each case. The ground state of L is dn(x,k), while the next two antiperiodic 
eigenfunctions are given by cn(x,k),sn(x,k). In what follows the parameter 
e = Vo + k 2 . There is an obvious rescaling which transforms the above form of 
the NLS equation to the form considered earlier. 

5.1 Focusing Case: 

It follows from Theorem (2) that for Vq + k 2 = e sufficntly small that the dn 
and sn solutions, as lower band edges, are modulationally unstable. Moreover, 
from Theorem (1) we are guaranteed a modulational instability as long as the L + 
operator is in the interior of a band at \x = 0. Thus in each case we have an open 
interval (0, e*) in which we are guaranteed a modulational instability, where e* is 
the smallest positive value of e such that fj, = is a band edge of the L+ operator. 
We will focus on the sn solution, since the instability of the dn solution follows 
from more elementary arguments []. We have found numerically that fi = is in 
the interior of a band for e e (-1.33, 0) U (-3.0,-1.56) U (-6.7, -5.6) . . ., which 
shows instability for e in these intervals. The numerical evidence further shows 
that in this particular case the second derivative of the branch of the Floquet 
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Figure 1: The Floquet Discriminants for the sn(x,k) solution in the focusing 
case with e = 0, .1 

discriminant passing through —2 is always negative, which implies that these 
solutions are always unstable. We do not currently have a proof of this. It is 
also unclear whether this feature is special to the elliptic function solutions, or if 
it holds in greater generality. It is interesting to note that the sn(x, k) solution 
also has a side-band type instability which appears at arbitrarily small positive 
amplitude. This is illustrated in|21 When e = the Floquet discriminants of the 
unperturbed problem fc(±/Lt) intersect near /i w .45. This is a collision of opposite 
Krein sign, since the discriminants have the same sign. This intersection opens 
into a gap under perturbation. 

The stability of the cn type solutions is extremely interesting. For small e 
the perturbation result guarantees that /i = is in a gap for the L + operator, 
and numerically we find no modulational instability for small e. As before we 
expect that for some critical value of e fi = will enter a band of L + and 
an instability will develop. In this case the critical value of e can actually be 
computed explicitly. When Vq = 0(e = k 2 ) we find that L+ is the 2-gap Lame 
operator, with fj, = is a band edge. It is also straightforward to show using 
the Sturm oscillation theorem and monotonicity of the band edges in e that 
this is the smallest value of e for which L + has a band edge at // = 0. Thus 
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we are guaranteed a modulational instability for some interval of beginning at 
Vq — 0(e = k 2 ). It is worth remarking that in this critical case (the integrable 
one) all of the information about the spectrum of the linearization, as well 
as a great deal more, can be obtained by the methods of algebraic geometry. 
See, in particular, Chapter 4 of the text by Belokolos, Bobenko, Enol'skii, Its 
and Matvccv 2 . The birth of a modulational instability at Vq = supports a 
intuition based on physical reasoning that was put forth in 3 which suggested 
that such solutions should go unstable at Vq = 0. 

Interestingly it appears that there exists a finite amplitude side-band type 
instability that sets in before the modulational instability. This is illustrated 
in figure [21 which shows the Floquet discriminants for the linearized operator 
of the focusing NLS about a cn type solution with modulus k 2 = 1/2 and a 
sequence of different values of e. Near Vq ph —.23 the Floquet discriminants 
K± cross, causing a gap to open and a loop of spectrum to emerge into the 
complex plane. This instability disproves a guess made in [3J, based on numerical 
experiments, that the cn type solutions should be stable for Vb < 0. As Vq 
increases the neighborhood about the origin in which the Floquet discriminants 
are real shrinks, until at Vq = the discriminants are only real at the origin. 
For positive Vo the discriminants again become real in a neighborhood of the 
origin. 

6 Conclusions 

We have established a sufficient condition for the modulational instability of 
a periodic standing wave solution to the NLS equation, which can be easily 
checked. In the case of weak nonlinearity this reduces to a physically reasonable 
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Figure 2: The Floquet discriminants for the focusing cn(a;, k) solution for Vq = 
-1/2, -1/4, -.228, -.1. 

criterion based on the effective mass of a particle in the periodic potential. 
We have also made some preliminary progress into studying the birth of side- 
band type instabilities. In the case of like Krein-sign collisions the gap-opening 
bifurcation is forbidden. Very preliminary numerical experiments suggest that, 
in the case of opposite Krein-sign collisions the "avoided collision" bifurcation 
is forbidden, but we do not currently have a proof of this. 
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an NSF FRG grant. 
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